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ABSTRACT 

Using our deep ~120 ks Chandra observation, we report on the results from our 
spatially-resolved X-ray spectral analysis of the "oxygen-rich" supernova remnant (SNR) 
0540-69.3 in the Large Magellanic Cloud. We conclusively establish the nonthermal na- 
ture of the "arcs" in the east and west boundaries of the SNR, which confirms the 
cosmic-ray electron acceleration in the supernova shock (B ~ 20-140 f^G). We report 
tentative evidence for Fe overabundance in the southern region close to the outer bound- 
ary of the SNR. While such a detection would be intriguing, the existence of Fe ejecta 
is not conclusive with the current data because of poor photon statistics and limited 
plasma models. If it is verified using deeper X-ray observations and improved plasma 



- 2 - 



models, the presence of Fe ejecta, which was produced in the core of the supernova, 
near the SNR's outer boundary would provide an intriguing opportunity to study the 
explosive nucleosynthesis and the ejecta mixing in this young core-collapse SNR. There 
is no evidence of X-ray counterparts for the optical O-rich ejecta in the central regions 
of the SNR. 

Subject headings: ISM: abundances — ISM: individual (SNR 0540-69.3) — supernova 
remnants — X-rays: ISM 



INTRODUCTION 



The supernova remnant (SNR) 0540-69.3 in the Large Magellanic Cloud (LMC) is well known 
for its bright central pulsar (PSR B0540-69.3) and its pulsar wind nebula (PWN) which show simi- 
lar characteristics to those of the Crab. The pulsar and PWN have been extensively studied at many 
wavelengths (e.g., Middleditch & Pennypacker 1985; Manchester et al. 1993a,b; Gotthelf & Wang 
2000; Kaaret et al. 2001; Petre et al. 2007 and references therein). SNR 0540-69.3 is also one of the 
two known members of "oxygen-rich" SNRs in the LMC. In the optical band, early data of 0540- 
69.3 showed a ring-li ke (~4 ;/ in radius) [O II I] feature around the pulsar/PWN, which is identified 



as the O-rich ejecta (jMathewson et al. 1980l ). The spatial and spectral structures of the optical fil- 
amentary features in this region were studied using the Hubble Space Telescope o bservations, which 
revealed the details of the ejecta and the synchrotron nebula (jMorse et al. 20061 ). The filamentary 
features of the optical line emission were suggested to result from magnetic Rayleigh- Taylor instabil- 
ities due to the interaction bet ween the expanding synchrotron nebula and the surrounding ejecta, 
just like those seen in the Crab (jMorse et al. 20061 ). This central region was also thorou ghly studied 
in infrared (IR) using the Spitzer Space Telescope observations (j Williams et al. 20081 ). These ob- 
servations revealed that the central nebula of 0540-69.3 is analogous to that of the Crab, showing 
synchrotron emission from the PWN, a complex network of ejecta filaments, and newly formed 
dust. A progenitor mass of 20-25 M was suggested based on the observed metal ejecta abun- 
dances ([Williams et al. 2 008). The enhanced optical and IR emission most likely originates from 
the interaction between the PWN and the freely-expanding metal-rich ejecta. 

Besides the central pulsar and PWN, the SNR shell produced by t he interaction between th e 
blast wave and the surrounding medium was detected in the radio band ((Manchester et al. 1993bl ) . 
The radio spectral index in the outer shell (a = 0.41, where S u tx v~ a ) is steeper than that in 
the PWN (a = 0.25), which indicated that the shell emission is synchrotron radiation from the 
shock-accelerated electrons. The age of SNR 0540-69.3 has been estimated to be t ~ 760-1660 yr 
based on various methods such as the pulsar spin-down, the kinematics of the optical ejecta, and 
the overall dynamics of the ejecta evolu tionary models (jSeward et al. 1984 : iKirshner et al. 19891 : 
Reynolds 19851 : iManchester et al. 1993al ). 



Thermal emission from the SNR shell (~1' in diameter) has been detected in the X-ray band 
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( Hwang et al. 200ll ; Ivan der Hevden et al. 200ll ). The overall SNR morphology of the nearly cir- 
cular shell with the bright W and faint E regions is consistent with that of the radio images. The 
bright W fil aments are roughly coi n cident with the patchy [O III] emission feature at ~30" W of 
the pulsar (jMathewson et al. 19801 ; ISeward k, Harnden 1994 ). The detection of these outer fila- 
ments, which most likely originate from interactions between the blast wave and ambient medium, 
in X-rays and radio makes 0540-69.3 observationally distinguished from the Crab in which no emis- 
sion features from the blast wave significantly beyond the PWN are detected at any wavelengths. 
Thermal emission from the outer SNR shell is detected prominently in X-rays, and thus X-ray 
observations of 0540-69.3 provide a unique opportunity to study the ambient and/or circumstellar 
structures in this SNR as well as the reverse shocked metal-rich ejecta. 

Previous X-ray studies of the SNR 0540-69.3 with Chandra and XMM-Newton produced puz- 
zling results. Enhanced line emission from O, Ne, and Fe was claimed by van der Hey den et al. 
(2001) based on the XMM -Newton RGS data, but was not detected based on the Chandra obser- 
vation (jHwang et al. 200ll ). The X-ray emission lines detected by XMM-Newton in the W regions 
of the SNR are strongly blue-shifte d (v ~ -2370 km s -1 , van der Heyden et al. 2001), which was 
not evident in the Chandra data (jHwang et al. 200ll ). The Chandra data suggested an intrigu- 
i ng possibility of th e presence of nonthermal X-ray emitting "arcs" in the E and W boundaries 



(jHwang et al. 200ll ). which are reminiscent of the nonthermal filaments of synchrotron radiation 
from the shock-accelerated cosmic-ray (CR) electrons in several Galactic SNRs such as SN 1006. 
However, the nature (thermal vs. nonthermal) of these faint arcs could not be determined with the 
previous Chandra data because of the poor photon statistics. The Chandra data also detected a 
hard tail (E > 3 k eV) in the observed X-ray spectru m which appears to prevail in the eastern hemi- 
sphere of the SNR (jHwang et al. 200ll : iHughes 20011 ) . The origin for the hard t ail (a thermal plasm a 
with kT ~ 5 keV vs. power law [PL] continuum with T ~ 2.7) was uncertain (jHwang et al. 200ll ). 
While the previous Chandra and XMM-Newton data revealed a number of new properties of SNR 
0540-69.3, some results were inconclusive and/or controversial. The uncertainties and discrepan- 
cies in the previous X-ray results appear to be caused largely by the poor photon statistics in the 
Chandra data (a 28 ks exposure) and/or the po or angular resolution (~1') of XMM-Newton RGS 



( Hwang et al. 200ll ; lvan der Hevden et al. 200ll ). Motivated by these discrepancies, we performed 
a deep Chandra observation to address the outstanding issues and to reveal the true nature of SNR 
0540-69.3. 

Here we report the results from the data analysis of our deep 120 ks Chandra observation of 
SNR 0540-69.3. We focus on our spectral analysis of the diffuse X-ray emission features outside of 
the PWN. The spectral analysis of the PWN will be presented elsewhere. In Section [21 we describe 
the observations and the data reduction. We present the spectral analysis of the SNR in Section [3j 
We discuss the results and implications in Section [H and present a summary in Setion [5j 
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OBSERVATIONS & DATA REDUCTION 



We observed SNR 0540-69.3 with the Advanced CCD Imaging Spectrometer (ACIS) on board 
Chandra X-Ray Observatory during A06 (in 2006 February 15-18, ObsIDs 5549, 7270, and 7271). 
We used the ACIS-S3 for the effective detection of the soft X-ray emission of the SNR. The previous 
Chandra observation of 0540 -69.3 (ObsID 119) showed significan t photon pileup effects in the 
central regions of the PWN (jHwang et al. 200ll : iPetre et al. 20071 ). To reduce the pileup effects 
in the PWN as much as possible while covering the entire SNR, we chose a 1/4-subarray of the 
ACIS. We corrected the spatial and spectral degradation of the ACIS data caused by the charge 
transfer inefficiency (CTI; Townsley et al. 2000) using the methods developed by Townsley et al. 
(2002). We applied the standard data screening by status and grade {ASCA grades 02346). Since 
our observation was split into three segments (~40 ks exposure for each segment), we performed 
these data reduction steps on individual observations. The three data sets were then merged into 
a single event file by re-projecting them onto the ObsID 5549's tangential plane. After the data 
reduction, the effective exposure was 114.2 ks. The archival Chandra data of 0540-69.3 (ObsID 
119) were obtained w ith a significantly shorter exposure (2 8 ks), and the results of the data analysis 
have been published (jHwang et al. 200ll : IPetre et al. 20071 ) . Since the low photon statistics of these 
archival data did not significantly enhance the statistics of our new data analysis, we used only our 
new deep observation in this work. 



3. Spectral Analysis 

Figure QJi shows an X-ray color image of SNR 0540-69.3. The central region is dominated by 
emission from the bright pulsar PSR B0540-69.3 and its PWN. The bright X-ray emission from the 
pulsar results in the trail image running NE-SW (corresponding to the CCD's read-out direction). 
We excluded this trail region (~2'/5 widt h) in the spectral a nalysis. SNR 0540-69.3 is projected on 



complex diffuse background structures (jHwang et al. 200ll ). We tested several nearby source- free 



regions for the background spectrum. Although regions to the S and W of the SNR show generally 
higher background emission (by ~20% in the total count rates) than regions to the N because 
of the relatively bright soft background emission (E < 1.5 keV) there, the selection of different 
background regions does not have a significant effect on our SNR spectral analysis. In this work, 
we use an average background spectrum extracted from four circular regions (~15" in radius) to 
the NE, NW, SE, and SW of the SNR (Figured]). 

Thanks to the use of the 1/4 subarray of the ACIS-S3, our new observation of SNR 0540-69.3 
shows reduced photon pileup effects in the PWN region. Based on simple tests using the event grade 
distribution and the count rate changes in the SNR between observations in 1999 and 2006, we find 
that regions with radial distance R ^ 3" from the pulsar are nearly free from pileup effects (Regions 
within R <^ 2" are still significantly piled up.) T he O-rich ejecta detected in the optical ban d shows 



a shell- like morphology with a radius of R ~ 4" (jMathewson et al. 19801 ; iMorse et al. 20061 ) . Thus, 
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using the pileup-free data, we can perform a reliable spectral analysis of this region to search for 
X-ray counterparts for the optical O-rich ejecta. However, the central regions of the SNR, which are 
spatially coincident with the optical ejecta features, are dominated by the bright PWN in X-rays. 
We find no evidence for enhanced X-ray emission features from the O He Q a nd /or Q Ly„ lines at R 
~ 4"— 5". This result confirms the conclusions inferred by the previous data (jPetre et al. 20071 ). We 
here concentrate our spectral analysis on the faint E-S regions, bright W filaments, and the faint 
E arc. In all spectral fits presented in the following sections, we binned each regional spectrum to 
contain a minimum of 20 counts per energy bin. 



3.1. The "Arc" Region 

Small arc-l ike filaments on the E and W boundaries of the SNR were discovered by the previous 



Chandra data ([Hwang et al. 20011 ). As Hwang et al. (2001) noted, this particular morphology 
and the hard X-ray spectrum of these features raised the intriguing possibility of strong particle 
acceleration sites at the shock front of 0540-69.3, in analogy to the "bi-polar" nonthermal X-ray 
filaments in the Galactic historical SNR 1006. In fact, Hwang et al. (2001) suggested that a 
significant portion of the X-ray emission in these arcs could be nonthermal continuum in origin. 
However, because of the poor photon statistics (e.g., ~190 counts in the arc region in the E 
boundary), the true nature of these features could not be revealed by the previous Chandra data. 
With a four times deeper exposure, we obtain significantly better photon statistics from these arcs 
to study the spectral characteristics. We note that the W arc region is projected on the bright 
thermal emission features of the SNR. Also, the trail image by the pulsar is partially superposed 
onto the W arc (Figure Q})) . The spectral analysis of the W arc is difficult and less useful than the 
E arc due to these complex environments and the detector artifacts. In fact, our two-component 
spectral model fits (plane shock + PL) for the W arc result in nearly the same spectral parameters 
for the hard PL component (r ~ 2.4, /0.5-8 keV ~ 4 x 10 -14 erg cm" 2 s _1 ) as those for the E arc 
(see below). Because the soft thermal emission dominates the overall W arc spectrum, the best-fit 
PL parameters are poorly constrained compared to those in the E arc. Thus, we discuss the nature 
of the arcs based on the spectral analysis of the E arc (the "Arc" region in Figure [2]) . 

The Arc region contains ~750 counts and is well-separated from the bright thermal filaments, 
as well as avoiding contamination by the pulsar trail image. The Arc region shows a featureless 
X-ray spectrum which can be fitted by a PL model with T = 2.42^q 2 3 (uncertainties are with a 
90% c onfidence level, hereafter) . We fixed the Galactic absorbing column at -/V^Gai = 7 x 10 20 



cm 2 (iDickey fc Lockman 19901 ). We fit the foregro und column in the LMC (A%,lmc = 6.5lj g x 



10 21 cm -2 ), assuming LMC interstellar abundances (IRussell Dopita 19921 ). The X-ray spectrum 
of the Arc region and the best-fit PL model (x 2 /v = 39.8/31) are presented in Figure [3J The 
plane-parallel shock model with LMC abundances cannot fit the observed spectrum (x 2 jv = 2.8). 
When the individual metal abundances are varied, the fit results in a large column (A^h,lmc ~ 
1.9 x 10 22 cm -2 ) and a high electron temperature (kT ~ 2.3 keV). The fitted abundances are 
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negligible for most species except for O (~5) and Ca (~4, abundances are with respect to Solar 
[Anders & Grevesse 1989], hereafter). The fitted LMC column is implausibly large for this location 
of the LMC. The inferred overabundances for O and Ca are statistically insignificant due to large 
uncertainties. The best-fit model (x 2 jv = 35.0/22) is significantly poorer than the PL model fit. 
A two-component model (PL + plane shock) indicates only a small contribution (~5% of the total 
flux) from the thermal plasma (kT ~ 0.4 keV). The PL + plane shock model does not improve 
the overall fit ( X 2 /v = 36.7/28). Thus, we conclude that the X-ray emission in the Arc region is 
dominated by nonthermal PL continuum. The best-fit photon index (r ~ 2.4) is consistent with 
the typical value for synchrotron radiation by relativistically accelerated electrons in the strong SN 
shock. 

Assuming that synchrotron radiation from the shock-accelerated electrons is responsible for the 
X-ray and radio emission for the Arc, we fit the observed Arc spectrum with the synchrotron cutoff 
model (SRCUT, Reynolds 1998; Reynolds & Keohane 1999). We fixed the radio spectral index a t 



a = 0.41 as estimated for the entire SNR excluding the central PWN ([Manchester et al. 1993bl ). 

hi.o 

-0.9 



The best-fit LMC column, A^hlmc = 5.7tn'o x 10 21 cm 2 , is consistent with that measured by 



the PL model fit. The fitted synchrotron roll-off frequency is v TO \\ = 1.01q 4 x 10 17 Hz. Our X-ray 
spectral analysis indicates that the 1 GHz flux for the radio counterpart is / iqhz ~ 0.19 mJy. On 



the other hand, based on the 1.5 and 5 GHz images ([Manchester et al. 1993b! ) . we roughly estimate 
/lGHz ~ 7 mJy for the Arc region. This radio flux is significantly larger than our estimate from 
the X-ray spectral fit, which suggests that the radio spectral index may not be uniform across the 
SNR. Thus, we repeated the SRCUT model fit of the Arc region by fixing /ighz = 7 mJy while 
varying the radio spectral index a. The best-fit parameters are -/V"h,lmc = 5.7j^g x 10 21 cm -2 , 
f ro ii = 2.lj^ J; x 10 17 Hz, and a = 0.62 ±0.01. The best-fit radio spectral index is not unusual for 



Galactic shell-type SNRs (jGreen 20091 ). Based on the PL and SRCUT model fits, the observed 
X-ray flux (absorbed) of the Arc is /o.s-8 keV ~ 3.7 x 10~ 14 erg cm -2 s . At the distance of d 
= 50 kpc for the LMC {d = 50 kpc for SNR 0540-69.3 is assumed in this paper hereafter), the 
unabsorbed X-ray luminosity is £o.5-8 keV ~ I- 7 x 10 er § s_1 - The best-fit model parameters for 
the Arc region are summarized in Table [TJ 

We note that the estimated absorbing column to the Arc can be reliably assumed for other 
regions of SNR 0540-69.3 because it is measured by straightforward models (PL and SRCUT) of the 
simple continuum-dominated spectrum. In fact, our spectral fits for the central region (R = 4"-5", 
the "ring" region in Figure [2]) also indicate -/Vh,lmc = 6.0±0.5 x 10 21 cm -2 . This central region 
is close to the PWN, in which the observed spectrum is dominated by the simple PL continuum 
and provides significant photon statistics (~3600 counts), while the pileup effects are convincingly 
ruled out. We find that other regional spectral analysis requires more complicated thermal spectral 
modelings which result in less reliable estimates of A^lmc (e-g-, Hwang et al. 2001). Thus, based 
on the results from the spectral analysis of the Arc region, we fix the LMC column at A^h.lmc 
= 6 x 10 21 cm" 2 in the following spectral analysis of the SNR. This A^h,lmc value is somewhat 
higher than those estimated/adopted in previous works (Nr ~ 4-5 x 10 21 cm -2 , e.g., Kaaret et al. 
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2001; Petre et al. 2007). The higher -/Vh,lmc is justified by the realistic assumption of the LMC 
abundances in our absorption model, while the previous works assumed the Solar abundances. Our 
-^h,lmc is fully consistent with that estimated by Serafimovich et al. (2004) who also assumed the 
LMC abundances for the absorbing column by the LMC. 



3.2. Faint East Region 

The eastern half of the SNR is generally faint, an d X-ray emission a ppears to originate from 



a hot plasma of the electron temperature kT ~ 5 keV (jHwang et al. 20011 ). However, the previous 
results on the faint eastern regions were rather preliminary based on a simple spectral analysis of 
the integrated spectra from the "inner" and "outer" regions of the entire eastern half of the SNR. 
The inner parts of the eastern half of the SNR could be substantially contaminated by scattered 
photons from the bright pulsar and PWN. The outer regions of the eastern half show an intensity 
variation with a relatively bright emission in the north. Furthermore, our new deep Chandra data 
reveal that X-ray emission near the southern boundary, which was included in the "outer east" 
region by Hwang et al. (2001), is spectrally-soft (Figure [1^) • Thus, we carefully re-analyze the 
faint E region excluding the relatively bright N regions and the spectrally-soft S regions. We also 
exclude the central regions of R ^ 20" to avoid the contamination from the pulsar and PWN as 
much as possible (our radial cut to remove the PWN is more conservative than that by Hwang et 
al. [2001] who excluded regions of R < 7".) Based on our MARX0 simulations of the pulsar /PWN 
of 0540-69.3, we estimated that the contamination by scattered photons from the pulsar/PWN due 
to the broad PSF wings of the mirror is ~20-30 % in the observed 3-7 keV band flux for the "E" 
region (as shown in Figure [2]). In these simulations, we assumed our observed ACIS image of the 
PW N (R < 4") with a PL spectrum of T = 2 and the 1-10 keV band flux of 3.2 x lO -11 erg cm" 2 



s _1 ( Hirayama et al. 20021 ). The spectrum of the scattered photons is hard (a PL with T ~ 0.25) 
at these large radial distances of R ~ 20"-30" from the pulsar. Thus, this instrumental artifact 
contaminates primarily at E ^ 5 keV. We fixed this scattered photon flux spectrum in the spectral 
model fits of the E region. (The broad PSF wing flux was estimated in the annular region of R = 
20"-30" using our MARX simulations, and was then normalized for the area of the E region. This 
normalized flux and the best-fit PL photon index for the scattered photon spectrum were fixed in 
the spectral model fits of the E region.) 

The X-ray spectrum of the E region includes ~1600 counts. The overall spectral characteristics 
are consistent with those shown for the "outer east" region in Hwang et al. (2001): the soft band 
spectrum evidently shows X-ray line emission features, and the hard tail extends to E > 4 keV. 
Initially, we fit the observed spectrum with a non-equilibrium ionization (NEI) plane-shock model 
(Borkowski et al. 200 1; vpshock in con junction with the NEI version 2.0 in the XSPEC) that is 



based on ATOMDB (jSmith et al. 200ll ). We use an augmented version of this atomic database 



^ttpV/space.mit.edu/CXC/MARX/ 
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to include inner-shell processes and updates of the Fe L-shell lines, whose effects are significant 
in highly under-ionized plasma spectrum, but are not accounted for in the current XSPEC NEI 
version 2.00 This model includes inner-shell Fe L lines from Na-like Fe ion at E = 725.15, 727.06, 
and 738.98 eV which are critical to properly fit the broad Fe L line complex-like feature in the 
0.7-0.8 keV band (see Section 3.3), but were not included in the standard XSPEC NEI version 
2.0. Thus, although it is still incomplete, our NEI model is the most reliable one that is currently 
available to the astronomical community to fit highly under-ionized plasma. We use this latest 
NEI model in the spectral analysis of thermal X-ra y emission throughout this work. We fixed 



the individual metal abundances at the LMC values (jRussell Sz Dopita 19921 ). We fixed the LMC 
column at iVj^LMC = 6.0 x 10 21 cm -2 based on the results from the spectral analysis of the Arc 
region (see Section f 3 . 1 [) . Although the overall spectrum can be fitted by a kT ~ 9 keV thermal 
plasma, the fit is statistically rejected (x 2 jv ~ 2.1) because of the large residuals in the soft band 
(E ~ 0.5-2 keV). When the LMC column is varied, the fit improves, but is still relatively poor 
(x 2 l v ~ 1-5). The fitted -/V^lmc (~ 2.5 x 10 21 cm~ 2 ) is significantly lower than that estimated 
for the Arc. Thus, the single plane-shock models with the LMC abundances are not appropriate 
to fit the observed spectrum. 

The fit improves to be statistically acceptable when the individual metal abundances are varied. 
We attempted to fit abundances for O, Ne, Mg, Si, S, Ar, Ca, and Fe while other elements (He, 
C, N, and Ni) were fixed at the LMC abundances. We found that only O and Ne abundances 
significantly affected the fit. Thus, we varied O and Ne abundances and fixed other elements at 
the LMC values (x 2 /^ = 64.8/55, Figure and Table [2]). The best-fit electron temperature is kT 
= 3.08l[j;92 keV with a low ionization timescale of n e t ~ 0.60 x 10 10 cm" 3 s. We fixed the LMC 
column at A^lmc = 6.0 x 10 21 cm -2 in this fit. The best-fit abundances for O (= 0.03^qq 2 ) and 
Ne (= 0.09^g Q3) are low, showing no evidence for metal-rich ejecta. When the LMC column was 
varied, the fit does not statistically improve ( x 2 /is = 62.8/54), and the best-fit JV H . LMC (= 1-2 X 
10 22 cm~ 2 ) is unreasonably high for the LMC. The best-fit abundances for O (~0.3) and Ne (~0.2) 
are similar to those of the LMC, but are poorly constrained. Thus, the E region spectrum can be 
described by a single plane shock model in which O and Ne abundances are lower than the LMC 
values. 

Alternatively, we fit the spectrum with a two-component plane shock model with the metal 
abundances fixed at the LMC values for both components. The LMC column was fixed at A^h,lmc = 
6.0 x 10 21 cm -2 for both components. The fit is statistically acceptable (x 2 /^ = 60.7/54, Table[2]). 
The best-fit electron temperatures are kT so { t = 1.72+Q3 3 and fclh ar d = 3.20t 2 'og keV. The soft 
component emission originates from an under-ionized thermal plasma (n e t S oft ~ 3.4 x 10 10 cm~ 3 s). 
The ionization timescale for the hard component is unconstrained. We also considered the hard tail 
of the observed X-ray spectrum to be nonthermal in origin, and thus repeated the two-component 



2 The augmented APEC atomic data have been provided by K. Borkowski. The relevant discussion on these 
detailed plasma model issues has been presented in Badenes et al. (2006). 
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model fits by replacing the hard component plane shock model with a PL. The best-fit PL photon 
index is V = 2.23±0.33, and the electron temperature of the thermal component is estimated to be 
kT = 1.72^069 (x 2 l v = 56.2/55, Table [2|). The modeled PL spectrum is consistent with that for 
typical synchrotron radiation from the shock-accelerated electrons. In fact, synchrotron radiation 
observed in the Arc region shows a similar PL spectrum with T = 2.4 (Section 13. ip . 

3.3. Soft South Region 

The faint southern boundary region, which is distinctively red in Figure is spectrally soft. 
This has not been recognized by previous data, and our new Chandra data reveal this feature. We 
extracted ~900 counts from this region (the "S" region in Figure [2]). The observed X-ray spectrum 
shows line emission features at E ~ 0.6-1.5 keV as well as the hard tail at E > 3 keV (Figuredb). 
We estimated the effects by the scattered photons from the pulsar and PWN due to the broad 
PSF wings of the mirror using the same methods described in Section 13.21 The contamination 
by scattered photons from the pulsar and PWN is estimated to be ~20% in the 3-7 keV band in 
this region. We fixed the PL component of the scattered photons (with the flux normalized for 
the area of the S region) in the spectral fits of the S region. Initially, we fit the spectrum with a 
single temperature NEI plane-shock model. We fixed the individual metal abundances at the LMC 
values, and also fixed the LMC column at iVf^LMC = 6.0 x 10 21 cm~ 2 . The best-fit model (kT ~ 
1.14 keV and n e t ~ 0.9 x 10 9 cm 3 s) is statistically poor (x 2 jv = 1.6). Varying iVn^MC does not 
improve the fit (-/V"h,lmc ~ 5.1 x 10 21 cm~ 2 , x 2 l v = 1-5). The poor fits are primarily caused by 
residuals at E ~ 0.7-0.9 keV. Because the residual line features at E ~ 0.7-0.9 keV appear to be 
the Fe L line complex, we varied the Fe abundance. The best-fit model (kT ~ 1.20^[J'gQ keV and 
n e t = 0.33±^;^ x 10 10 cm" 3 s) statistically improves (x 2 jv = 32.9/33, Figuregb and Tabled]), and 
suggests overabundant Fe (= 18.5^3 j). We similarly attempted to vary O and/or Ne abundances, 
and found that O and Ne overabundances are not required to improve the fit. When O and/or Ne 
abundances are fitted with the Fe abundance fixed at the LMC value, the best-fit model fit clearly 
shows residuals at E ~ 0.7-0.9 keV, which suggests that the emission features in the 0.7-0.9 keV 
band are likely dominated by the broad Fe L line complex. Varying A?h,lmc results in nearly the 
same spectral parameters with no statistical improvement (x 2 jv = 32.3/32): The Fe overabundance 
is still suggested (Fe = 17.3 j^) with similar kT (~1.37 keV) and 7V h ,lmc (~5.6 x 10 21 cm" 2 ). 

For completeness, we attempted two-temperature plane-shock models to fit the S region spec- 
trum. ]Vh,lmc was fixed as above, and the metal abundances were also fixed at the LMC values. 
The best-fit model (x 2 '/v = 31.4/31) suggests a high temperature plasma (^TWa > 2.4 keV) in 
addition to the soft component (kT so { t ~ 0.35 keV). Because the single temperature model fit 
suggested an Fe overabundance, we tested whether the Fe overabundance is required in the two 
component model. When the Fe abundance is varied, an Fe overabundance (>1.2) is suggested, 
while other spectral parameters remain generally consistent with those estimated with the LMC 
abundances (x 2 /v = 25.8/30). Nonetheless, these two-component models do not statistically im- 
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prove the fit, compared to the single component models. Based on the single plane shock model 
fits, the enhanced Fe abundance appears to required to adequately fit the observed spectrum of 
the S region. However, the statistical uncertainties on the Fe abundance measurements are large 
(i.e., the Fe overabundance is a ~2-3cr detection above the LMC value). Also, the two-temperature 
plasma models may equally fit the observed spectrum without overabundant Fe. Thus, while the Fe 
overabundance is suggested in the S region, a firm conclusion may not be drawn with the current 
data. 



3.4. Bright West Regions 

In contrast to the faint E and S regions, the western parts of the SNR show bright fila- 
mentary emission features (Figure [Th.). Hwang et al. (2001) revealed that the bright western 
filaments, except for the small region near the outer boundary in the direction opposite to the E 
arc, are dominated by soft X-ray emission from thermal plasma with kT ~ 0.5-1 keV. The previous 
Chandra data indic ated no evidence for enhanced metal abundances in the bright western regions 



(jHwang et al. 20011). while the XMM-Newto n RGS data suggested enhanced emission from Fe L 



and Ne K lines (|van der Heyden et al. 20011 ). Results from our spectral analysis of the western 
filaments are generally consistent with those by Hwang et al. (2001). We here present our spectral 
analysis of two representative regions: the "SW" and "NW " regions (Figure [2j). We choose the 
SW region because it shows bright thermal emission at the outer boundary, likely representing 
the blast wave with good photon statistics. The NW region includes a significant fraction of the 
hard emission (kT ~ 5 keV) among the western filaments (besides the regions of the candidate 
nonthermal filament at the outer boundary, Hwang et al. 2001). In the spectra of the SW and 
NW regions, the contamination by the scattered hard X-ray photons from the pulsar and PWN is 
negligible thanks to the small extraction area. Thus, we do not include the PL component of the 
broad PSF wings in the spectral fits of these regions. 

The SW region spectrum (~2200 counts) is presented in Figure Eh. The overall spectrum can 
be fitted by a single component plane-shock model (kT ~ 0.74 keV, n e t ~ 5.3 x 10 10 cm -3 s, x 2 /v 
= 74.5/54). We fixed the LMC column at A^h^mc = 6.0 x 10 21 cm -2 , and the individual metal 
abundances at the LMC values. When A^h,lmc was fitted, the best-fit value (./Vh,lmc ~ 6.9 x 
10 21 cm -2 ) is similar to that estimated from the Arc, and the overall fit does not improve (\ 2 jv = 
74.0/53). Varying the individual abundances results in the LMC-like abundances (~0. 1-0.3) with- 
out a significant improvement in the fit (x 2 jv = 60.8/49) based on an F-test (F-probability ~ 0.07). 
There is no evidence for the presence of a hard component (E ^ 3 keV) or metal overabundances 
in this region. The best-fit spectral parameters assuming the LMC abundances are summarized in 
Table [2J 

For the NW region spectrum (~2800 counts), a single plane-shock model fit with the fixed 
-^h,lmc and the LMC abundances is not acceptable (x 2 jv ~ 2.0). When the LMC column was 
fitted, a ~50% higher column (Ah.lmc = 9.1^0 x 10 21 cm" 2 , kT ~ 1 keV, \ 2 /v ~ 1.6) is implied. 
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However, the observed spectrum of this region is nearly identical to that of SW region in the soft 
energy band {E <^ 1 keV, Figure H]), in which the large foreground absorption in the NW region, if it 
exists, should significantly affect the observed spectral shape. Then, we varied the metal abundances 
(with -/Vj^lmc = 6.0 x 10 21 cm -2 ), which results in a significantly improved fit (x 2 /v = 80.1/65, 
Figure Hti). The best-fit abundances are low (~0.05-0.2) showing no evidence for metal-rich ejecta. 
The fitted spectral parameters are summarized in Table [2j Since t he hard spectral co mponent 
with kT ~ 5 keV was suggested in this region in the previous work (jHwang et al. 200ll ). we, for 
completeness, attempted a two-component plane-shock model fit for the NW region. In this fit, 
we fixed A^lmc = 6.0 x 10 21 cm -2 and the abundances at the LMC values. The best-fit model 
indicates kT so f t ~ 0.6 keV and kT^d ~ 1-4 keV (x 2 / u ~ 1-0)- Although the two temperature shock 
model can adequately fit the observed spectrum, the electron temperature for the hard component 
is significantly lower than that implied in the E region. Varying A^h,lmc an d the abundances does 
not affect these results. The best-fit spectral parameters for the two component model are also 
summarized in Table [2j 



4. Discussion 



4.1. Nonthermal Filament 



Using our deep Chandra observation, we firmly establish the nonthermal nature of the Arc in 
the E boundary of SNR 0540-69.3. The observed X-ray spectrum is well fitted by a PL with T 
= 2.42, which is typical for synchrotron radiation from relativistically accelerated electrons in the 
strong shock. The contribution from underlying thermal emission is estimated to be small (~5% of 
the total flux, if exists). Assuming the synchrotron origin, our SRCUT model fits indicate that the 
roll-off frequency between the radio and X-ray bands is ~ 1-2 x 10 17 Hz. The Chandra image 
shows that the Arc consists of several substructures with various widths and brightness. Based on 
these filamentary features, we estimate the widths of the Arc to be in a range of ~l"-2 / . / 8 (FWHM) 
which correspond to a physical size of ~0.3-0.7 pc. SNR 0540-69.3 may be in a free-expansion 
phase at the age of r ~ 1000 yr (e.g., Hwang et al. [2001] and references therein). Then, the 
angular distance of the Arc to the SNR center (~35" corresponding to a physical distance of ~8.5 
pc) suggests a time-average shock velocity v s ~ 8300 r^ 1 km s -1 , where t\ is the SNR age in units 
of 1000 yr. Assuming the Sedov phase, v s ~ 3300 t^qq km s" 1 is implied. These velocities are 
in plausible agreement with the estimated shock velocities (v s ~ 2000-9000 km s _1 ) for the SNR 
in a free-expansion or a Sedov phase, based on the spectral analysis of the thermal component 
emission for the E region (Section 14. 2p . The estimated widths of the Arc filaments and these shock 
velocities (t> s ~ 2000-9000 km s _1 for the E region) imply the synchrotron loss time ti oss ~ u v^ 1 r 
~ 160-2700 yr, where u is the advection distance of the downstream electrons or the physical width 
of the filament, and r is the compression ratio in the shock. In this estimate, we assumed r ~ 5-8 
for an efficient particle acceleration (e.g., Ellison et al. 2007). 
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The peak frequency of a synchrotron emitting electron is u p = 1.8 x 10 18 E% B± Hz, where 
B± is the magnetic field component perpendicular to the velocity vector of the electron, and E e is 
the electron energy. For v p ~ 3 x 10 17 Hz (or the photon energy ~1.3 keV) representing typical 
X-ray photons as observed in the Arc spectrum, the corresponding electron energy E e = 0.4 B±~2 

_ 2 

erg implies B± = 136 ti oss 3 G. For the estimated ri oss ~ 160-2700 yr, we derive B± ~ 7-46 jiG. 
Considering a geometrical projection effect in estimating the widths of the filaments, the actual 
B field is likely larger than these estimates. For instance, the observed filament width would be 
overestimate d by a factor of ~5 assuming a spherical shock with an exponential emission profile 
(jBallet 2 006). Then, B ~ 20-140 /iG are estimated. The maximum electron energy is E e>max = 
2.5 x 10 -7 f ro ii 5 B^g 1 TeV ~ 7-25 TeV. The estimated B field is roughly an order of magnitude 
higher than typical LMC interstellar fields (B ~ 4 /xG, Gaensler et al. 2005), indicating a strong 
magnetic field amplification in the SNR shock. It was suggested that the presence of ordered 
magnetic fields in the LMC (and in other young/irregular galaxies) could have originated from the 
CR-driven d ynamo due to recent SN activity rather than the standard galactic rotational dynamo 
mechanism (jGaensler et al. 2 005). The detection of synchrotron emission by the shock- accelerated 
relativistic electrons in SNR 0540-69.3, which is located in the LMC's largest active star-forming 
region 30 Doradus, provides support for the presence of such a population of relativistic particles 
in the LMC. These spectral characteristics of the nonthermal emission, the morphology of the 
thin filamentary features in opposite sides of the SNR, and the young age of ~1000 yr for SNR 
0540-69.3 ar e similar to thos e observed in Galactic SNRs 1006 (e.g., Bamba et al. 2003) and 
G330. 2+1.0 (|Park et al. 20091 ) . While this general scenario is plausible, we note that the quality of 
the current data is somewhat limited, and thus our simplified interpretations of the magnetic field 
amplification need to be taken with a caution. For instance, the detailed spatial substructures of the 
Arc is complex rather than a spherical shell. Thus, our simple geometrical correction may not be 
fully justified, and the true uncertainties on the magnetic field measurements could be significantly 
larger. High angular resolution X-ray imaging data with deeper observations will be required to 
conclusively reveal the true nature of the magnetic field structure in this SNR. 



4.2. Thermal Emission and Metal Abundances 

Based on the spectral analysis of the E region, we confirm the presence of a hard component 
(E > 3 keV) in 0540-69.3. Two-component spectral model fits indicate that this hard component 
emission can be described either by a hot gas {kT ~ 3 keV) or a PL (r ~ 2.2). With the current 
data, it is difficult to statistically discriminate the thermal or nonthermal origin for the hard X-ray 
emission in the E region. Although the PL component is similar to that for synchrotron radiation 
in the Arc region (r ~ 2.4), there is no evidence for thin filamentary features which would have 
been anticipated along the outer boundary of the E region, had the hard X-ray emission originated 
from the shock accelerated relativistic electrons in highly amplified magnetic fields just behind the 
shock front. However, we note that SNR 0540-69.3 is significantly more distant than the Galactic 
SNRs that show thin nonthermal filaments. If the E region of 0540-69.3 consisted of a nest of 



-13- 



thin filaments as observed in the Galactic nonthermal SNR G347.3-0.5, it could be seen largely 
as a faint smooth rim except for the brightest filaments in the Arc region. Thus, the nonthermal 
interpretations for the hard component emission in the E region are not clear with the current data. 
On the other hand, the hard tail of the E region spectrum can be described by a hot gas of kT 
~ 3 keV (either by a single or two component shock model fits). The presence of this hot plasma 
is plausible in young SNRs. Although a nonthermal origin for the hard X-ray emission in the E 
region cannot be ruled out with the current data, we here, for simplicity, discuss the nature of the 
E regional spectrum assuming a thermal origin. 

Based on the best-fit volume emission measure (EM = n e n R V = 1.74 x 10 57 cm" 3 ), estimated 
by the simple one-component plane-shock model fit, we derive the post-shock electron density n e 
~ 1.44 / 2 cm , where / is the volume filling factor of the X-ray emission. In this estimate, we 
assumed n e = 1.2nn (where «h is the H number density). We also assumed the X-ray-emitting 
volume V ~ 1.0 x 10 57 cm 3 for the E region with an average path-length of ~2.4 pc along the 
line of sight, corresponding to the angular thickness of the E region (~10") at the distance of d 

_ 1 _„ 

= 50 kpc. The two-temperature model implies n e ~ 0.67 and 1.38 / 2 cm for the soft and the 
hard component, respectively. These estimated post-shock n e are several times lower than those 
measured in the bright W regions (see below). The estimated ionization timescales (n e t ~ 0.6-3 
x 10 10 cm~ 3 s) indicate the ionization ages t ~ 130-1600 yr depending on models. The swept-up 
mass in this region, M sw = nwin p V (where m p is the proton mass), is estimated to be small (M sw 
~ 1.0-1.4 /3 Mq), which is consistent with the shock propagating through a low-density medium 
in a young SNR. 

In contrast, the X-ray spectrum in the bright W regions can be fitted by thermal plasmas with 
relatively low electron temperatures (kT ~ 0.5-1.5 keV). These spectra also imply low LMC-like 
metal abundances. Thus, the bulk of X-ray emission in the W portions of the SNR originates 
from the shocked interstellar /circumstellar medium. Based on the fitted EM, we estimate the 
post-shock n e for the SW and NW regions. In these calculations, we adopt the same method as 
above. For the SW region, we use the X-ray emitting volume V ~ 9 x 10 55 cm -3 , assuming the 
path-length of ~1 pc along the line of sight. The estimated post-shock electron density is n e ~ 
5.6 /~2 cm~ 3 . For the NW region, we assume V ~ 5 x 10 55 cm~ 3 for the path-length of ~0.8 pc 
along the line of sight. The estimated density is n e ~ 6.8-12.3 /~2 cm -3 depending on models. 
The blast wave is encountering several times denser ambient medium in western regions than in 
eastern regions. While the surrounding density is significantly different between E and W sides of 
the SNR, the overall radial extent is nearly symmetric in E-W as well as in N-S. We speculate that 
the ambient density is higher in W, and that the bla st wave has prob ably entered the dense medium 



there recently. Although the HI map of the LMC (jKim et al. 19981 ) shows that SNR 0540-69.3 is 



in fact surrounded by relatively bright HI filaments, these HI filaments are primarily located in N 



and E sides of the SNR. The bright W shell is not detected in the mid-IR (IWilliams et al. 2 008). 



On the other hand, Hwang et al. (2001) noted that the CO map of the LMC (j Cohen et al. 19881 ) 



shows molecular clouds to W of the SNR. However, the low-resolution of the CO map (~8'8) does 
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not allow one to determine whether the CO clouds are close enough to the SNR to account for 
the higher implied density in the western filaments. Alternatively, the dense material in the W 
regions of the SNR might be a relic structure of the stellar winds from the massive progenitor. 
In this scenario, however, the stellar winds must be highly asymmetric between E and W, which 
is difficult to reconcile with the winds from a 20-25 M & progenitor in the red supergiant phase. 
High-sensitivity, high-resolution radio observations of 0540-69.3 regions would be useful to study 
the structures of the surrounding medium and its interactions with the SNR, which may help reveal 
the origin of the bright W shell of the SNR. 

If the SNR is in the adiabatic phase, we can derive the Sedov parameters of the SNR us- 
ing the results of our spectral analysis. We attempt these estimates using the characteristic E 
region with a representative electron temperature of kT = 3 keV. Assuming an electron-ion tem- 
perature equilibrium, this temperature implies a shock velocity v s ~ 1600 km s _1 . However, the 
assumed electron-ion te mperature equilibrium is unlikely established for v s ^ several 10 2 km s _1 



( Ghavamian et al. 20071 ). Thus, ~1600 km s _1 may only be considered as a lower limit. The radial 
distance from the pulsar to the outer boundary in the E region is ~30" corresponding to ~7.3 pc, 
and thus the Sedov age of r se d < 1800 yr are estimated for v s > 1600 km s _1 . This SNR age limit 
is in good agreement with the previously estim ated SNR ages of ~ 760-1660 yr. Alternatively, the 



SNR might still be in the free-expansion phase ([Hwang et al. 200ll ). especially in the eastern parts 



where the density is low. Assuming the free-expansion phase, the shock velocity should be v s ~ 
7150 rf 1 km b -1 . If we take tsnr ~ 800 and 1800 yr as lower and upper limits for the SNR age 
appropriate for the spectral parameters derived in the E region, the corresponding shock velocities 
are v s ~ 4000-9000 and 1600-3600 km s" 1 for a free -expansion and a Sedov phase, respectively. 
On the other hand, in the SW region, the electron temperature of kT = 0.74 keV implies v s > 800 
km s _1 . The radial distance from the SNR center to the SW boundary is also ~30", and thus r se d 
< 3600 yr is estimated. While this age range is broadly consistent with the previous estimates, 
the derived upper limit on the SNR age is larger than that estimated in the E region. Considering 
the fact that the angular extent is similar between the E and W boundaries, the large Sedov age 
estimated in western regions may support that the blast wave has recently entered a dense medium 
in western regions as discussed above. 

The spectrally-soft emission in the S region appears to be caused by enhanced emission from 
the Fe L line complex (the one-shock model fits) and/or the low temperature of the plasma (the 
two-shock model fits). A low electron temperature (kT ~ 0.35 keV) for this region is implied by 
the two-component shock model fits. This temperature is ~5-10 times lower than those estimated 
in the E region, which could suggest a significantly higher density and thus substantial brightening 
in the X-ray intensity in the S region compared with its surroundings, unless the low temperature 
component originates from small clumpy material that may not be resolved by the current data. 
However, the overall X-ray emission in the S region is faint with no evidence for large density and 
high X-ray intensity, which is generally not consistent with the inferred low temperature. In the 
two-shock model fits, the existence of the hard component (kT > 2.4 keV) is not firmly conclusive 
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due to the limited photon statistics (only ~40 source counts at E > 3 keV). In fact, the one- and 
two-shock model fits are statistically indistinguishable (Section I3.3[) . and the presence of the two 
characteristic plasma components may not be convincingly justified in the current data. 

The single shock model fit suggests Fe overabundance in the S region. Although the detection 
of the Fe overabundance in the S region is not conclusively established in the current data, if true, 
the suggested candidate Fe-rich ejecta near the southern outer boundary of the core-collapse SNR 
is intriguing. The presence of the Fe overabundance in the S region would imply that the Fe-rich 
ejecta created at the deepest core of the massive progenitor has been expelled out to the outermost 
boundary of the SNR. Such a case has been detected in the young Galactic O-rich SNR Cas A 
(e.g., Hughes et al. 2000, Hwang & Laming 2003). Although this general scenario is intriguing, a 
firm detection of the Fe ejecta in the S region is required before discussing conclusive and extensive 
interpretations of this candidate Fe ejecta. High photon statistics with deeper X-ray observations 
and improved NEI models to accurately account for the significantly under-ionized Fe L lines are 
required to perform a quantitative abundance study of this soft feature in the S region, and thus 
to reveal the true nature of this candidate Fe ejecta. 

We do not find evidence for enhanced line emission from Ne and Fe in the bright W regions. 
This is inconsistent with the results reported b y van der Heyden et al. (2001) using the XMM- 
Newton RGS data (jvan der Heyden et al. 200 ll ). and the discrepancy remains unanswered. Van 
der Heyden et al. (2001) claimed moderate overabundances (by a factor of ~2-3) in Ne and Fe 
relative to the O abundance in the W regions, without quoting uncertainties. We note that our 
Chandra data hint for less abundant O than Ne and Fe (e.g., in the bright NW region, we estimate 
the O abundance is ~20% of the LMC value, while the Fe abundance is ~50% of the LMC). Thus, 
we speculate that high abundances for Fe and Ne relative to O based on the XMM-Newton data 
could have been due to low-abundant O rather than the enhanced Ne and Fe. Finally, we find no 
evidence for the O-rich ejecta corresponding to the optical ring of [O III] enhancements close to 
the PWN. For the annular region of R = 4"-5" (the "ring" region in Figure [2]), the observed X-ray 
spectrum can be fitted by a PL model with T = 2.39±0.09 (A^lmc = 6.0±0.5 x 10 21 cm" 2 , xl 
= 163.3/123). We place a 3a upper limit of / < 5.2 x 10" 16 erg cm" 2 s" 1 (= 5.0 x 10" 7 photons 
cm" 2 s" 1 ) for the O Ly a line flux. 



5. Summary 

We performed a deep ~120 ks Chandra observation of the O-rich SNR 0540-69.3 in the LMC. 
We conclusively establish the nonthermal nature of the "arc" in the E boundary of the SNR. The 
observed X-ray spectrum shows a PL continuum (r ~ 2.4), indicating synchrotron radiation from 
the shock-accelerated relativistic electrons. Thus, CR electron acceleration is evidently effective 
in the shock front of this young SNR. Assuming that synchrotron radiation is responsible for the 
observed X-ray and radio emission in this nonthermal filament, we estimate a roll-off frequency of 
z-Wi ~ 1~2 x 10 17 Hz. From our estimates of the shock velocities (v 8 ~ 2000-9000 km s" 1 ), the 
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physical width of the filament, and the observed spectral shape, we derive post-shock magnetic 
fields B ~ 20-140 /j,G. The maximum electron energy is estimated to be £ e ,max ~ 7-25 TeV. 

We confirm the presence of the hard component emission in this young SNR primarily in the 
faint eastern regions. The hard component spectrum can be described either by a hot gas with kT 
~ 3 keV or a PL with T ~ 2.2. It is difficult to discriminate the origin of the hard component 
(thermal vs. nonthermal) because of the low photon statistics in the current data. If the hard 
X-ray emission is nonthermal, the estimated PL photon index suggests synchrotron radiation for 
its origin, which is similar to that detected in the Arc region. If the hard emission is thermal, 
the hot gas implies an ambient density of no ~ 0.17-0.35 cm -2 which is several times lower than 
that estimated in the western regions (no ~ 1.4-3 cm -2 ). The W regions are dominated by bright 
soft filamentary emission features representing the shocked high density medium of lower electron 
temperatures (kT ~ 0.5-1.5 keV). Metal abundances in the W regions are low (^ 0.3), indicating 
interstellar and/or circumstellar origins rather than metal-rich stellar ejecta. Our spectral analysis 
of thermal X-ray emission features indicates high shock velocities (v 8 ~ 2000-9000 km s -1 ) which 
are usually expected in young SNRs where particle acceleration is efficient. 

In the S boundary of the SNR, we detect a faint emission feature which shows a distinctively 
soft X-ray spectrum compared to the surrounding regions. The Fe abundance in this region appears 
to be enhanced, although the detection of the Fe overabundance is not conclusive because of the low 
photon statistics in the currrent data and the incomplete Fe L line models for highly under-ionized 
thermal plasma. If it is confirmed by follow-up high statistics data and improved NEI models, the 
suggested Fe-rich ejecta candidate at the outer boundary of the SNR is intriguing. It would be the 
detection of Fe ejecta created in the deepest core of the SN and subsequently expelled to nearly 
the outer blast wave just like those seen in Cas A. Deeper Chandra observations and improved Fe 
L line models for significanly under-ionized plasma are required to reveal the true nature of this 
candidate Fe-rich ejecta and its implications on the core-collapse nucleosynthesis and the details of 
the explosion mechanism. 
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Table 1. Best-Fit Model Parameters for the Arc Region. 



Model 


Ah.lmc 
(10 21 cm" 2 ) 


^roll 

r (io 17 Hz) 


1 GHz Flux 
(mJy) 


a 


x 2 /^ 


PL 
SRCUT 
SRCUT 


6-5±i;| 

5-7^ 
5-^3 


9 a 9+0.26 
z -^ z -0.23 

1 +L4 

2 1+ 4 ' 7 

z - i -1.2 


0.19±0.02 
7 b 


0.41 a 
0.62±0.01 


39.8/31 
41.6/31 
41.1/31 



a The radio spectral index a (wher e S v tx v a ) is fixed at a = 0.41 as estimated for 
the entire SNR excluding the PWN ([Manchester et al. 1993bl ). 

b The 1 GHz flux is fixed at /igHz = 7 mJy as we estimate based on the 1.5 GHz 
image by Manchester et al. (1993). 



Table 2. Best-Fit Spectral Parameters of SNR 0540-69.3. 



Region 


feTaoft 

(keV) 


&7hard 

(keV) 


r 


^cisoft 

(10 10 cm" 3 s) 


^eihard 

(10 10 cm" 3 s) 


EMsott 
(10 57 cm" 3 ) 


-EMhard 

(10 57 cm" 3 ) 


x> 


E a 

S b 
SW 
NW C 


3.08l°;f 2 

1 79+1-63 
l-' z -0.35 
i 79+2.3S 
l--' -0.69 

1 2fl +0 19 

1 - zu -0.60 
u -'*-0.09 

qn +0 16 

u - d9 -0.04 


o 9fi+2.55 
°- zu -1.06 

-, 4 o+0.75 


2.23+0.33 


60+ 012 

Q O7+3.03 

— 1.47 

3 4fi+ 616 
^ +0 ' n 

u -°°-0.07 
r 97+3.43 
<->- z ' -1.67 

6.30!£» 
50.7tf- 6 


<0.08 


1 74+0-47 
!■ '^-0.27 

0.37l°;J 3 
0.30t°;i» 

1 14+0-40 
1.1*_0.23 
9 Q9+0-35 
Z - OZ -0.26 
n qo+3-96 
u.oo_ 1 gg 

4 00+0-48 
"±.»^_0.63 


1 5Q+0.52 

1 qq + 1-70 
J--3J-0.85 


64.8/55 
60.7/54 
56.2/55 
32.9/33 
74.5/54 
80.1/65 
66.3/68 



Note. — -/Vh.lmc is fixed at 6.0 x 10 21 cm 2 . Errors are with 90% confidence. For each region, one- and two-component 
plane-shock model fits are presented. The metal abundances are fixed at the LMC values, unless noted otherwise. 

a The O and Ne abundances are varied in the one-shock model. 

b The Fe abundance is varied while other elemental abundances are fixed at the LMC values. 

c The metal abundances are varied in the one-shock model, while abundances are fixed at the LMC values in the two-shock 
model. 
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Fig. 1. — (a) An X-ray color image of SNR 0540-69.3 taken from the 114 ks Chandra observation. 
Red is 0.4-0.7 keV, green is 0.7-3 keV, and blue is 3-7 keV bands. Each subband image is binned 
into 0"492 pixels, (b) The hard band (2.5-7 keV) Chandra image of SNR 0540-69.3. The image is 
binned into 0'.'492 pixels. For the purposes of display, the image is smoothed by a Gaussian with a 
= 1". In both panels, the central pulsar and PWN are saturated to white to emphasize the faint 
filamentary features of the SNR. 
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Fig. 2. — A grey-scale broadband (0.4-7.0 keV) Chandra image of SNR 0540-69.3. The image 
binned into / .'492 pixels. Regions of our spectral analysis are marked. Dashed circles show regions 
where the background spectrum is estimated. 
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Fig. 3. — The Chandra ACIS spectrum of the Arc region. The best-fit PL model is overlaid. The 
lower panel shows the residuals from the best-fit model. 
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Fig. 4. — The Chandra ACIS spectrum of (a) the E, (b) the S, (c) the SW, and (d) the NW regions. 
In (a)-(d), the thick solid line shows the best-fit single plane-shock model. The lower panel shows 
the residuals from the best-fit model. In (a) and (b), the dotted line shows a PL model representing 
the estimated contamination by the scattered photons from the PWN. 



